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li.  '  ABSTRACT 


5  A  description  is  given  of  the  development  and  corroboration  of  a 
simplified  computational  model  for  the  prediction  of  the  radiated  rotational 
noise  of  a  lifting  rotor  or  propeller.  The  method  is  based  on  a  solution  of 
the  concentrated  force-excited  wave  equation  and  the  identification  of  the 
terms  in  this  solution  with  annular  distributions  of  monopole  sources  of 
specified  phase  and  amplitude. 

The  computational  algorithms  developed  from  this  mathematical  model 
provide  a  rapid  means  for  determining  the  amplitude  and  phase  of  the  rad¬ 
iated  sound  field.  They  are  particularly  well  suited  for  providing  a  des¬ 
cription  of  rotor  noise  characteristics,  which  can  be  used  as  input  to  com¬ 
puter  programs  designed  to  calculate  the  rotor  noise  field  in  the  presence 
of  boundaries.  (  ' 

Results  obtainedwith  this  method  usually  differ  by  less  than  4  dB 
from  results  obtained  by  the  analytical  method  of  Lowson  and  Ollerhead. 
Comparisons  of  the  theory  with  the  experimental  cruise  data  of  Schlegal,  et 
al,  and  with  results  of  whirl-tower  tests  by  Stuckey  and  Goddard  is  generally 
satisfactory.  In  those  cases  where  the  predicted  and  measured  data  disagree 
by  more  than  5  dB,  it  is  felt  that  this  difference  can  be  explained  by  the 
lack  of  accurate  blade-loading  data  and  by  experimental  error. 
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Preface 


This  report  summarizes  the  results  of  a  study  to  develop  a 
mathematical  model  and  computational  procedure  for  character¬ 
izing  the  rotational  ..-use  of  a  helicopter  rotor.  The  resulting  data 
would  be  suitable  as  input  for  computer  programs  designed  to 
calculate  the  rotor  noise  field  in  the  presence  of  boundaries.  This 
work  was  performed  by  the  Ordnance  Research  Laboratory  of  The 
Pennsylvania  State  University  under  contract  IsOw  65-0123-d  for 
the  Aeronautics  Group,  Code  461,  of  the  Office  of  Naval  Research, 
Arlington,  Virginia. 


Abstract 


DESCRIPTION  is  giver,  of  tne  development  and  corroboration 
of  a  simplified  computational  model  for  the  prediction  of  the  radi¬ 
ated  rotational  noise  of  a  lifting  rotor  or  propeller.  The  method 
is  based  on  a  solution  of  the  concentrated  force-excited  wave 
equation  and  the  identification  of  the  terms  in  this  solution  with 
annular  distributions  of  monopole  sources  of  specified  phase  and 
amplitude. 

The  computational  algorithms  developed  from  this  mathemat¬ 
ical  model  provide  a  rapid  means  for  determining  the  amplitude 
and  phase  of  the  radiated  sound  field.  They  are  particularly  well 
suited  for  providing  a  descrij  tion  of  rotor  noise  characteristics, 
which  can  be  used  as  input  to  computer  programs  designed  to 
calculate  the  rotor  noise  field  in  the  presence  of  boundaries. 

Results  obtained  with  this  method  usually  differ  by  less  than 
4  dB  from  results  obtained  by  the  analytical  method  of  Lowson 
and  Ollerhead.  Comparisons  of  the  theory  with  the  experimental 
cruise  data  of  Schlegal,  et  al,  and  with  results  of  whirl-tower 
tests  by  Stuckey  and  Goddard  is  generally  satisfactory.  In  those 
cases  where  the  predicted  and  measured  data  disagree  by  more 
than  5  dB,  it  is  felt  that  this  difference  can  be  explained  by  the 
lack  of  accurate  blade-loading  data  and  by  experimental  error. 


Nome*  clatu*  e 


Symbol 


.(r.0) 


A 


escriptii  i 


Coefficient  of  Fourier  series 
describing  blade  spatial 
pressure  variation! 

Complex  pressure  amplitude 
of  monopole  source  at  unit 
distance 


Am<r-®> 

B 

c 

c0 

F(r,©.t) 

«  » 

•  • 

t  F'(r.e.t) 

l  " 

|  i. 

? 

i  — 


Unit-distance  amplitude  of 
source  monopole  at  r,9  for 
mth  sound  harmonic 

Number  of  blades 

Chord  length 

Speed  of  sound 

Complex  amplitude  of  point 
force  expressed  in  rotor 
polar  coordinates  for 
impulsive  chordwise  loading 

Complex  amplitude  of  force 
expressed  in  rotor  polar 
coordinates  for  rectangular 
chordwise  loading 

Retarded  force 
Radial  force  component 
Rotor  (hub)  altitude 


-  <^gsns -«W>*.  -aw*SS5»i«SSK®K> 


h  (r,@) 
m ' 


h,n(x-y'  z> 


J 

k 

k. 


*1  «^2>^3 

L(r.0) 


M 


n(or  \) 

N. 


p(x,  y,  z.t) 


Pme(x.y.z) 


Fr 


•  rO 


q(t) 

q#(t) 


Coefficient  of  Fourier  series 
for  the  blade  force  expressed 
in  rotor  Dolar  coordinates 


Coefficient  of  Fourier  series 
for  the  blade  force  expressed 
in  rectangular  coordinates 


.1/2 


(-D 

Effective  ring  factor 


■ju 

c0 


Direction  cosines  of  the 
force  vector  coefficient 


Blade  section  loading 


Sound  harmonic  number 


Hub  Mach  number.  — 

cn 


Air-load  harmonic  number 


Number  of  dipoles  in  acoustic 
array  model 


Total  acoustic  pressure  at 
(x.y,  z.t) 


Acoustic  pressure  of  mth 
line  component  at  (x.y.  z) 


Dipole  pressure  at  "r"  (see 
Fig.  3) 


Pressure  at  "r"  due  to  monopole 
at  point  "0"  (see  Fig.  3) 


Pressure  at  ’’r"  due  to  monopole 
at  point  "1"  (see  Fig.  3) 


Mass  flux  rate  of  monopole 


Time  derivative  of  mass  flux 
rate  of  monopolc 


Rotor  torque 


Rotor  polar  coordinate  (s-c 
Fig.  1) 


i  ! 


B  ? 


u„.  Vr 


x,  y,  z 


ryr  'i 


x..  y.,  z. 

i  i  i 


X..  Y.,  Z. 

ill 


x  ,  y  ,  z 
r  r  r 


Rotor  or  propeller  radius 


Instantaneous  distance  between 
moving  monopole  and  receiver 


[(yr-ys)2+(2r-zs)2l1/2 
(see  Fig.  5) 


[(x-xi)2+(y-y1)^+(z-z1)2]1/2 
(see  Fig.  1) 


Sound  pressure  level 
(abbreviation) 


Time  variable 


Rotor  thrust 


Normalized  Fourier  coefficients 
of  blade  load  series 


Hub  velocity 


General  field  point  coordinates 


X  coordinate  of  moving  monopole 
at  retarded  time 


General  source  coordinates 


Coordinates  of  monopole 
member  of  ith  dipole 


Coordinates  of  monopole 
member  of  ith  dipole 
corrected  for  rotor  tilt 
angle 


X  coordinate  of  moving 
monopole  at  time  "t" 


Receiver  coordinates 


Blade  angle  of  attack 


Angle  between  source  velocity 
vector  and  vector  joining 
the  source  and  receiver  at 
retarded  time  for  a  moving 
mcnopole 


An  unspecified  phase  angle 


Spacing  of  monvpoles  for 
implementing  dipole  radiation 
characteristics 


‘•^**^<5 *>&!tt?%<‘VUf  *  \Syt*  f^.'  ; 


A0 

6 

®, 

X  (or  n) 
’n 

$c 

* 


Dipole  angular  separation 

Azimuthal  angle  (see  Fig.  1) 

Mean  azimuth  angle  for  ith 
dipole 

Air-load  harmonic  number 

Phase  angle  for  nth  air-load 
harmonic 

Blade  coning  angle 

Rotor  tilt  angle 

Source  radian  frequency 


Rotor  angular  velocity 
(rpm  x  2t?/60) 


n  mBS 

‘  m 


1 


An  Acoustic-Array  Model  for  the  Computation 
of  the  Rotational  Noise  of  a  Lifting  Rotor 

The  increasing  importance  of  civilian  and  military  VTOL  craft 
(especially  helicopters)  has  revived  interest  in  the  estimation  and 
reduction  of  propeller-  or  rotor-generated  noise;  however,  the 
only  rigorous  theory  that  exists  applies  to  the  noise  from  a  pro¬ 
peller  w,th  symmetrical  sections,  zero  blade  angle,  and  zero 
forward  speed  (Ref.  1).  There  are,  however,  various  approximate 
analyses  that  can  be  used  to  estimate  rotational  noise  for  practical 
thrust- producing  rotors  and  propellers. 

The  earliest  and  still  most  commonly  employed  approach  to 
estimating  propeller  noise  is  based  on  a  method  originally  em¬ 
ployed  by  Gutin  (see  Ref.  2).  Basically,  this  method  consists  of 
postulating  an  approximate  pressure  distribution  over  the  blades 
and  then  using  this  distribution  to  calculate  the  characteristics 
of  independent  fixed  forces  that,  when  distributed  over  the  area 
swept  by  the  blades,  produce  an  acoustically  equivalent  radiation 
pattern.  The  radiated  sound  field  is  computed  as  the  solution  to 
the  appropriate  wave  equation,  assuming  the  energy  source  is 
adequately  represented  by  the  fixed-force  distribution. 

Various  extensions  and  refinements  of  Gutin’s  method  hive 
been  reported  in  the  literature  (see,  for  example.  Refs.  3,  4,  and 
5).  These  refinements  appear  to  provide  the  basis  whereby  rea¬ 
sonably  accurate  estimates  of  rotor  noise  may  be  obtained,  at 
least  for  the  free-field  rotational  noise  of  isolated  propellers  or 
helicopter  rotors  and  providing,  of  course,  that  proper  consider¬ 
ation  is  given  to  the  effects  of  nonuniform  blade  loads. 

A  more  recent  method'  for  computing  rotor  rotational  noise 
stems  from  the  work  of  Lowson  (Ref.  6)  who,  instead  of  basing 
his  development  on  the  concept  of  fixed  forces  distributed  in  the 
blade  path,  approaches  the  problem  by  considering  the  sound  fields 
arising  from  singularities  in  motion.  Specifically,  equations  are 
derived  that  describe  the  general  acoustic  relationships  for  a 
point  force  in  arbitrary  motion  as  well  as  those  for  a  point  source 
(acoustic  monopole)  and  point  acoustic  stress  in  arbitrary  motion. 
Moreover,  it  is  shown  how  these  relationships  can  be  applied  to 
the  calculation  of  various  aspects  of  aircraft  noise.  These  results 
reduce  to  those  of  Gutin  for  the  case  of  a  uniformly  loaded  pro¬ 
peller.  i 


In  piratical  situations,  it  is  often  unsatisfactory  to  regard 
propellers  or  rotors  as  isolated  noise  sources  due  to  the  baffling 
effects  of  the  vehicle  fuselage  or  appendages.  Of  even  greater 
importance,  in  some  cases,  is  the  necessity  of  obtaining  the  noise 
field  of  VTOL  craft  not  in  the  free  field  but  rather  in  the  vicinity 
of  an  extended  boundary  such  as  the  ground  or  ar.  ocean.  In  view 
of  these  considerations,  this  study  was.  undertaken  to  develop  a 
mathematical/ computational  model  of  a  helicopter  rotor  (pro¬ 
peller)— a  model  that: 

1.  Can  be  more  readily  applied  to  the  estimation  of  rotational 
noise  fields  in  the  presence  of  boundary  surfaces 

2.  Can  provide  a  simplified  method  for  computing  the  free-field 
rotational  noise  of  an  isolated  rotor  for  comparison  with 
other  data. 


* 


i 


Development  of  the  Mathematical  Model 


* 

1 


fig.  1  -  Rotor  geometry. 


The  metnud  developed  here  results  in  a  description  cf  the 
rotor-system  rotational  noise  in  terms  of  an  equivalent  array  of 
independent,  acoustic-monopole  (point),  time-harmonic  sources. 
Evidently,  this  characterization  is  considered  complete  once  the 
spatial  location,  amplitude,  and  phase  of  each  source  in  the  array 
are  specified. 

The  principal  assumptions  made  in  the  course  of  the  mathe¬ 
matical  analysis  are: 

1.  Blade  angle-of-attack  is  independent  of  azimuth 

2.  The  total  torque  and  thrust  of  the  blade  system  are  considered 
to  be  concentrated  within  a  thin  annular  ring  located  in  the 
path  of  the  blades  and  at  a  radius  equal  to  some  fraction  k 
of  the  blad":  radius 

3.  The  thrust,  torque,  and  radial  force  vectors  operate  in  pnase 
with  one  another  at  any  particular  point  in  the  blade  plane 

4.  The  blade  chordwise  pressure  distribution  is  assumed  to  be 
rectangular  and,  in  the  limiting  case  of  an  infinitely  narrow 
blade,  impulsive 

5.  Compressibility  effects  are  negligible 

6.  Transient  effects  are  ignored 

7.  The  acoustic  medium  is  unbounded,  homogeneous,  and  iso¬ 
tropic.  The  effect  of  boundaries  is  handled  separately  by  con¬ 
sidering  the  reflection  and  transmission  of  sound  at  existing 
boundaries  for  each  simple  source  in  turn  and  then  combining 
these  fields  by  using  the  superposition  principle. 
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Fig.  2  -  Pressure-vs-tirce  history  it 
tr.  e.  h*rO,t. 


SPATIAL  AND  TEMPORAL  PRESSURE 
VARIATIONS  IN  ROTOR  DISK 

Figure  1  shows  the  rotor  geometry  used  in  the  analysis.  The 
rotor- shaft  axis  is  represented  by  the  y  axis;  the  time  and  blade 
azimuthal  references  are  referred  to  the  x  axis. 

Figure  2  shows  an  assumed  rectangular  pressure-«ersus- 
time  history  experienced  by  a  stationary  observer  located  at  a 


i 


S 


i 


\ 

k 

t 

2 

■5 

"M 

3 

l 

z 


■i--  .-'y2-2j:  >  »>:-A  Zi^r?** 


point  in  the  path  of  the  rotor.  From  this  figure  it  follows  that  the 
pressure  can  be  written  as  a  Fourier  series: 


Note  that  the  m  =  0  term  has  been  neglected,  since  it  does  not 
produce  an  acoustic  signal.  In  terms  of  blade-section-loading 
data,  i.e.. 


L(r,0)  =  cP(r,e)  , 


we  have 


...  ,  ®  2L(r,@)  .  /m8c\  „/  0  C  \ 

(2)  Pir.9.n.z  s“>(ir)cosmflB(,'n'1flr)  ■ 


The  magnitude  of  the  force  corresponding  to  this  pressure  is 
obtained  by  multiplying  Eq.  2  by  r-dr-d  0: 


(3) 


F'(r.0,t) 


*  2rL(r,0) 

L  — — —  Sin 

m®l 


mile 


(f) 


cos  mas(t  -  -  ~)  drd©  . 


The  rectangular  pressure  distribution  assumed  in  the  deriva¬ 
tion  of  the  preceding  equation  leads  to  results  that  tend  to  deem- 
phasize  or  even  eliminate  certain  rotor-noise  line  components 
due  to  the  presence  of  the  factor  (l/m)sin(mBc/2r).  This  behavior 
is  not  generally  observed  in  experimental  data;  however,  precise 
definition  of  the  chordwise  pressure  distributions  can  be  expected 
to  be  significant  only  at  frequencies  at  which  the  wavelength  is 
not  greater  than  the  order  of  a  chord  length  (see  Ref.  7).  Accord¬ 
ingly,  for  most  practical  cases  that  arise  in  rotor-noise  calcula¬ 
tions,  one  finds  that  a  point-loading  distribution  (impulse)  is 
adequate  for  the  purpose  of  representing  the  actual  distributed 
chordwise  pressure  distribution.  In  addition  to  its  mathematical 
simplicity,  the  impulsive  loading  assumption  has  the  advantage 
that  it  leads  to  results  that  appear  to  be  more  nearly  reflected 
by  available  experimental  data.  Other  chordwise  pressure  dis¬ 
tributions  can  be  accounted  for  by  reevaluating  Eq.  1. 

An  expression  equivalent  to  Eq.  3  for  point  (impulsive)  loading 
can  be  easily  derived  in  the  limiting  case  of  small  chord  length. 
In  fact,  since 


(1) 

P (r,@,  t)  =  £  om  (r,8)  cos  mftBt', 
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«  • 
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- 
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lim  — - - —  =1  , 
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it  follows  from  Eq.  3  that  the  expression  describing  the  amplitude 
and  phase  of  the  point  force  at  ( r,  8,  t )  is 

(4)  F(r,0,t}=  £  L{r,  6)  cos  mflB^t- drd@  . 

ACOUSTIC  RELATIONS 

We  begin  with  the  force-driven  wave  equation,  which  relates 
the  acoustic  pressure  to  the  driving-force  vector: 


151 


In  the  limit  of  a  point-applied  force,  the  solution  to  Eq.  5  is  (see 
Ref.  8) 


S’Mx-x,)*  +  (  y-y,)2  +(z-z,)z  . 


When  the  force  vector  at  a  point  in  the  rotor  disk  can  be  rep¬ 
resented  as  a  Fourier  series,  such  as 


(7)  F(r.©,0  =  I  hm(r,®)  cos(ilm1  +  y)  , 

m»l 


then  the  solution  (Eq.  6)  may  be  written  as 

(8,  P(,.,.vi-4 1 

Now,  since  we  can  write  the  force-vector  coefficient  in  terms  of 
its  magnitude  hm(x,y,z)  and  its  direction  cosines  f,,f2,and  f3, 
Eq.  8  becomes 
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4 


and  y  is  an  arbitrary  initial  phase  angle.  From  a  comparison  of 
Eqs.  10  and  11  it  is  evident  that  the  pressure  field  in  Eq.  10  is 
made  up  of  monopole  sources  for  which 
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A»(r,©) 
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We  can  now  relate  Eq.  12  to  Eq.  4,  which  was  derived  to  de¬ 
scribe  the  amplitude  of  the  force  whose  line  of  action  is  along 
the  vector  8S.  For  a  force  directed  to  p.-oduce  a  net  lift  on  a 
horf-ontal  rotor  system,  as  in  Fig.  1,  we  find  that 

(13)  h*(r,8)  =  L  (r,8)  drd@  , 

and  the  line  of  action  of  the  force  vector,  as  well  as  SS,  can  be 
taken  to  be  normal  to  the  blade  chord.  If  we  now  make  use  of  the 
assumption  that  the  rotor  forces  are  concentrated  within  an 
“effective  ring”  located  at  a  radius  equal  to  kR,  we  can  avoid 
the  radial  integration  implied  by  Eq.  13  and,  using  a  harmonic 
series  for  the  normalized  blade-loading  data,  write 


tw  L(,,w*.[i'.(^),.r:]  (ttMM 


cos  [n(0-8r 


It  is  noteworthy  that  the  effective- ring  approximation  has  long 
been  used  in  the  application  of  propeller-ncise  theory,  and  it  has 
the  effect  of  simplifying  the  computations —  usually  without  intro¬ 
ducing  severe  errors.  With  regard  to  its  use  in  predicting  heli¬ 
copter  noise,  Lowson  and  Ollerhead  (Ref.  7)  point  out  that,  under 
the  assumption  of  random  loading-phase  variations  around  the 
rotor  azimuth  and  along  the  blade  span,  the  error  arising  from 
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consideration  of  only  one  radial  station  (i.e.,  a  single  effective 
ring)  is  small,  at  least  in  the  far  field.  A  representative  figure 
for  the  location  of  a  single  effective  ring  is  taken  to  be  the  80 
percent  radial  station  (k=0.8). 

The  amplitude  of  the  monopole  sources  comprising  the  equiv¬ 
alent  rotor  noist  model  of  the  mth  line  component  of  rotational 
noise  can  be  computed  by  combining  Eqs.  12,  13,  and  14: 

I/I 

(15)  |A.(kR.0)|«-^-{l+l[uJ+v:] 

411*8 S  L  "*|L  J 

cos[n(©-©t  )~CB]  }  A©  , 

where 

1/2 

F'  =  [T'*(iSf)'*F’]  • 

and  where  the  symbol  for  the  infinitesimal  d©  has  been  replaced 
by  the  source  angular  separation  A@.  Thus,  it  is  proposed  to 
model  the  acoustic  field  of  the  lifting  rotor  by  means  of  two 
circular  arrays  of  phase-  and  amplitude-shaded  monopoles  ar¬ 
ranged  as  shown  in  Fig.  4.  Here,  the  upper  ring  of  sources  is 
composed  of  positive  monopoles  whose  amplitude  and  phase  are 
related  to  the  azimuthal  angle  assigned  to  each  monopole  by  the 
expression 

F.  r  «  T  l'/I  r  n  Fig.  4  -  Uampleo;  m!ci  acoustic-array  model. 

{,6)  A-(kR’0)=^i?H, luH  Hnl0-0'KJ} 

(a  \ 

,-lT/  0  ■ 

The  lower  monopole  ring  is  composed  of  sources  of  opposite 
sign.  Note  that  the  two  source  rings  are  slightly  displaced  from 
one  another  in  azimuth  (and  have  slightly  difiarent  diameters)  so 
that  the  dipole  vectors  are  aligned  with  the  blade  force  vectors. 

The  acoustic  field  of  the  mth  line  component  developed  at  a 
field  point  with  coordinates  x,  y,  z  due  to  the  operation  of  the 
monopole  at  x^,  y.,  (kit,  ©  in  rotor  polar  coordinates)  follows 
immediately  from  Eq.  16: 

■/> 

cos  . 

In  Eq.  17  the  plus  sign  is  used  for  monopoles  in  the  upper  ring 
and  the  minus  sign  for  those  in  the  lower  ring.  The  total  acoustic 


pressure  due  to  all  the  sources  in  the  array  is  obtained  simply 
by  superimposing  the  effects  of  the  individual  sources  at  the  re¬ 
ceiver.  This  computation  is  readily  adapted  for  machine  imple¬ 
mentation. 


SPATIAL  RELATIONS 

The  basic  rotational  noise  model  has  been  described  in  terms 
of  an  array  of  phase-  and  amplitude-shaded  acoustic  monopoles, 
it  remains  to  establish  the  number  and  spatial  location  of  the 
sources  in  the  array. 

It  is  convenient  to  specify  the  location  of  the  source  rings  so 
that  the  mean  rotor  height  (h)  is  midway  between  the  two  y  coordi¬ 
nates  describing  the  upper  and  lower  ring  heights.  If  the  azimuthal 
angle  8  is  referred  to  the  horizontal  plane  located  at  the  mean 
rotor  height  (and  with  8=0  in  the  direction  of  the  positive  x  axis), 
the  location  of  the  monopoles  may  be  specified  using  the  blade 
angle  of  attack  a  t  the  coning  angle  <pc,  and  the  mean  azimuth 
angle  8|  corresponding  to  the  ith  dipole.  In  fact,  the  coordinates 
of  the  two  monopcles  comprising  the  ith  dipole  are: 

z  SS  r  i 

(18)  x,  =r  cosdit-^-^s'masinB.-cos  cc  cos8,  sin  <f>c  j 

t  t,SS 

y ,  =  ht— cos  a  cos  <#>e 

zf- -r  sin8,t-^  [since  cos  8,+ cos  a  sin  8,  sin  <#>c  ] 

where  the  upper  sign  is  used  for  the  sources  in  the  upper  ring 
and  the  lower  sign  for  the  sources  in  the  lower  ring. 

The  case  of  a  rotor  (or  propeller)  arbitrarily  oriented  with 
respect  to  the  assumed  coordinate  system  is  handled  by  straight¬ 
forward  coordinate  transformation  once  the  equivalent  acoustic 
array  is  computed  for  the  horizontal  rotor  configuration.  The 
transformation  employed  here  uses  a  rotor  tilt  angle  which 
specifies  the  inclination  of  the  rotor  disk  relative  to  the  x-z  plane. 
The  positive  x  axis  is  used  as  the  angle  reference  and  an  upward 
tilt  of  the  rotor  plane  is  defined  as  positive.  Parenthetically,  it 
should  be  mentioned  that  it  is  considered  unnecessary  to  provide 
for  rotation  of  the  blade  system  about  the  rotor  axis,  since  the 
same  relative  effect  can  be  obtained  by  selecting  appropriate 
field  point  (receiver)  coordinates.  Thus,  the  transformation 
equations  may  be  written  as: 

(I9>  X;=(xf)cos  ^-[(y7.)-h]sin  + 


Yf  =  h  +  (x7)sm  +  [(y^)-h  ]  cos  * 
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To  retain  adequate  fidelity  in  the  process  of  computationally 
simulating  the  rotor  noise  sources,  it  is  necessary  to  provide  a 
sufficient  number  of  dipoles  in  the  approximation  to  the  rotor 
noise  continuum.  The  basic  consideration  is  one  of  using  enough 
dipoles,  equally  spaced  around  the  effective  ring,  to  satisfy 
Shannon’s  sampling  theorem  for  the  highest  air-load  harmonic 
to  be  considered. 

The  problem  of  estimating  the  range  of  air-load  harmonics 
required  for  the  computation  of  the  rotational  noise  of  the  mth 
line  component  was  considered  by  Lowson  and  Ollerhead  (Rrf.  7). 
They  concluded  that,  for  accuracy,  air-load  data  must  be  included 
in  the  range  given  by  the  expression 

(20)  mB(M-l)<  X  <  mB(M  +  l) 

Accordingly,  to  provide  an  adequate  number  of  dipoies  in  the 
simulation  of  the  highest  air  load,  it  is  necessary  that 

(21)  Ns>2mB(M+l) 


For  the  third  acoustic  line  component  of  the  SH-3A  (HSS-2) 
helicopter  rotor,  for  example,  Eq.  21  specifies  that  about  44 
dipoles  (88  monopoles)  should  be  used  for  the  highest  required 
air-load  harmonic  (the  22nd). 

Another  quantity  that  must  be  provided  when  using  the  model 
is  the  magnitude  of  the  monopole  spacing  $S .  Generally,  it  has 
been  found  that  the  relation 


(22) 


8S  =0.01 


2Tko 
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gives  good  results  computationally. 


The  Computational  Model 

The  computational  model  consists  of  algorithms  for  specifying 
the  location  (xj,  yj,  zj  coordinates),  amplitude  (at  unit  distance), 
and  relative  phase  of  each  source  in  an  acoustic  array  designed 
to  simulate  the  rotor  noise.  Equations  16  and  19  are  the  basic 
expressions  that  describe  these  quantities. 

From  Eq.  16,  which  gives  the  strength  and  phase  of  a  monopole 
in  the  circular  array,  we  note  that  rather  detailed  knowledge  of 
the  blade  loading  distribution  is  assumed.  Unfortunately,  blade 
loading  characteristics  are  not  well  known  at  present.  For  ex¬ 
ample,  no  information  is  available  that  will  permit  one  to  calculate 
the  phase  angle  .  Furthermore,  information  on  the  blade-load 
amplitudes  (Un  and  Vn)  is  empirical  and  is  based  chiefly  on  the 
in-flight  measurements  reported  by  Scheiman  (Ref.  10),  which 
were  conducted  using  an  instrumented  helicopter. 

Using  Scheiman's  data,  together  with  that  of  Burpo  and  Lynn 
(Ref.  11),  Lowson  and  Ollerhead  (Ref.  7)  proposed  a  loading  law 
for  the  blade-load  amplitudes  based  on  the  inverse  2. t>  power  of  the 
blade- loading  harmonic  number  X  assuming  use  of  the  effective- 
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ring  approximation.  To  circumvent  the  specification  of  individual 
phase  angles  they  further  assumed  that  the  acoustic  contribu¬ 
tions  of  the  various  air-load  harmonics  could  be  combined  in  the 
mean-squared  sense. 

The  proposed  i  coustic-array  geometry  that  emerges  as  a 
result  of  the  above  considerations  is  composed  of  a  number  of 
circular  arrays  similar  to  that  shown  in  Fig.  4.  Due  to  the  fact 
that  the  £„  are  as  yet  unknown,  however,  a  single  equivalent 
array  for  the  rotor  noise  field  cannot  be  obtained.  Instead,  it  is 
necessary  to  employ  individual  circular  arrays  to  account  for 
the  contributions  from  each  of  the  air-load  harmonics  considered 
and  to  sum  the  resultant  acoustic  fields  in  the  mean-squared  sense 
as  suggested  above. 
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CORRECTIONS  FOR  ROTOR-SYSTEM  TRANSLATION 
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Fig.  5  -  Source-netion  geometry. 


It  is  convenient  to  account  for  the  effects  of  rotor-system 
translation  by  applying  the  proper  uniform-motion  corrections 
to  the  pressure  field  of  each  monopole  before  superimposing 
their  individual  contributions  at  the  receiver.  For  the  case  of  a 
stationary  observer  located  at  coordinates  (xr,  yr,  zr),  and  for 
uniform  subsonic  hub  motion  in  the  direction  of  the  Dositive  x 
axis,  the  appropriate  corrections  can  be  derived  beginning  with 
the  results  of  Morse  and  Ingard  (Ref.  12).  After  a  few  typograph¬ 
ical  errors  in  Eq.  11.2.15  of  the  above  reference  are  corrected, 
one  can  write  this  equation  as 

I  q[t-~^-}cos0-M)V 

(23)  p  - - +  - 

4n  R,(I-M  cos  j91*  4IIRrI(l-Mcos/9)5 

This  equat  >n  expresses  the  observed  pressure  from  a  moving 
inonopole  in  terms  of  the  rate  of  mass  flux  out  of  the  source 
(i.e.,  its  strength  q(t),  the  time  derivative  of  the  mass  flux 
rate  (q'  (t) ),  the  source-to-receiver  distance  at  retarded  time 
(Rr),  and  the  angle  between  the  source  velocity  vector  and  the 
vector  joining  the  source  and  receiver  at  retarded  time  ( /3)  . 

For  a  harmonic  monopole  source  whose  pressure  field  at  unit 
distance  is  given  by  the  expression 

P=A  e1"’  , 
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Eq.  23  takes  the  form 

(24)  p  =  — - - ; 

Rr(l-Mcos/3) 


(cosjS-M)M  { 
*  *  k0  R,(l-M  cos/3)  J 


To  compute  the  appropriate  values  of  Rr  and  (3,  we  make  use 
of  the  geometry  shown  in  Fig.  5.  We  note  that,  without  loss  of 
generality,  the  motion  of  the  source  may  be  taken  to  be  along  a 
line  parallel  to  the  x  axis  and  at  a  distance  zs  from  the  y  axis. 
The  instantaneous  position  of  the  source  is  given  by  the  coordi¬ 
nates  (xs,  ys,  zs);  t"e  source  position  when  the  sound  was  emitted 


Thus,  given  the  location  (xs,  ys,  zs),  as  well  as  the  amplitude  and 
phase  (A)  of  the  source,  Eq.  24  can  be  used  to  compute  the  con¬ 
tribution  of  the  source  at  the  stationary  receiver  position  for  any 
subsonic  value  of  the  huo  convection  Mach  number,  M=Y/cg. 


Illustrative  Examples 


The  first  test  of  the  mathematical  model  was  made  by  com¬ 
paring  various  computational  results  with  those  obtained  using 
the  method  suggested  by  Lowson  and  Ollerhead  (Ref.  7).  To  main¬ 
tain  reasonably  close  correspondence  between  the  two  techniques, 
the  present  computations  were  also  carried  out  using  the  assump¬ 
tions  embodied  in: 

1.  A  loading  law  for  the  harmonic  amplitudes 

2.  A  10:1:1  relationship  between  the  magnitudes  of  the  force 
vectors  for  T,  Q,  and  Fr 

3.  A  random  phase  relationship  betw-een  the  acoustic  contribu¬ 
tions  of  the  various  air  loads 

4.  The  effective- ring  approximation. 

Note  that  these  assumptions  are  not  indigenous  to  the  present 
computational  procedure.  In  fact,  the  first  two  assumptions  can 
be  eliminated  by  the  simple  expedient  of  providing  the  computer 
program  with  data  cards  that  describe  some  other  set  of  loading 
amplitudes,  phases,  and  force-component  magnitudes;  whereas 
the  remaining  two  can  be  eliminated  by  making  only  minor  changes 
in  the  computational  algorithms. 

Figure  6  is  a  polar  plot  of  the  predicted  sound  patterns  for 
the  first  three  line  components  of  a  Sikorsky  S-58  helicopter 
rotor  in  the  hover  configuration.  These  sound  patterns  were 


Fig.  6  -  Free-lield  sound  patterns  ol  the 
Sikorsky  S-5S IHSS-H  helicopter  in  hover 
configurelion. 
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computed  for  a  tadial  distance  from  the  hub  of  320  ft;  therefore, 
they  are  approximately  representative  of  free-field  conditions. 
The  figure  also  includes  the  number  of  azimuthal  stations  (number 
of  dipoles)  as  well  as  the  range  of  air-load  harmonics  (X)  used 
by  the  digital  program.  The  very  close  agreement  between  the 
results  of  the  two  techniques  is  evident. 

Figure  7  shows  the  predicted  sound  patterns  for  four  line  com¬ 
ponents  of  the  SH-3A  rotor  in  the  hover  configuration.  Again  the 
data  are  in  close  correspondence.  The  predicted  sound  patterns 
for  four  line  components  of  the  same  rotor  at  a  radial  distance 
much  nearer  to  the  hub  and  essentially  in  the  n»ar  field  of  the 
rotor  is  shown  in  Fig.  8.  At  this  distance  the  far-field  assumptions 
inherent  in  the  method  of  Lowson  ar.i  Ollerhead  are  probably 
unjustified;  nevertheless,  it  can  be  seen  that  the  data  produced 
by  use  of  the  two  methods  are  in  close  agreement. 


Fitj.  7  -  Free-field  sound  patterns  ot  the 
Sikorsky  SH-3A  (HSS-21  helicopter  in  hover 
configuration. 


Fig.  8  ■  Near-field  sound  patterns  of  the 
Sikorsky  SH  3A  (HSS-21  helicopter  in  hover 
configuration. 
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a  In  general,  use  of  the  present  computational  procedure  to 

accurately  predict  the  near-field  noise  of  a  helicopter  rotor 
requires  that  reliable  data  on  the  loading  phase  variations  be 
used  and  that  the  assumption  of  an  effective  ring  be  discarded. 
Nevertheless,  it  appears  that  useful  data  can  be  obtained  in 
many  cases  even  if  these  refinements  are  not  made.  In  particular, 
this  conclusion  seems  to  be  valid  for  field  points  outside  the 
t  „  sphere  defined  by  the  path  of  the  rotor  tip. 

Figures  9,  10,  and  11  illustrate  the  calculated  effects  of 
forward  speed  on  the  sound  patterns;  they  show  the  instantaneous 
sound  patterns  for  the  SH-3A  rotor  at  steady  forward  flight  speeds 
of  120,  80,  and  40  knots.  For  comparison,  several  representative 
data  points  calculated  by  the  method  of  Lowson  and  Ollerhead  are 
provided  for  the  120-knot  case  (Fig.  9).  Again  the  two  methods 
. .  yield  similar  results. 
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fig.  9  -  filed  of  forward  speed  on  the  sound 
patterns  of  the  SH-3A  helicopter  rotor  in  steady 
flight  at  a  forward  speed  of  120  knots. 
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Fig.  11  -  Effect  of  forward  speed  on  the  sound 
patterns  of  the  SH-3A  helicopter  rotor  in  steady 
flight  at  a  forward  speed  of  40  knots. 
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Fig.  10  ■  Effect  of  forward  speed  on  the  sound 
patterns  of  the  SH-3A  helicopter  rotor  in  steady 
flight  at  a  forward  speed  of  80  knots. 
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To  test  the  present  theory  against  measured  data,  sound 
pressure  levels  coi  responding  to  the  “flyby”  (cruise)  measure¬ 
ments  reported  by  Schlegal,  et  al,  were  calculated  for  the  S-58 
helicopter.  These  data  correspond  to  the  40-,  80-,  and  110-knot 
cruise  cases  and  for  the  first  four  acoustic  line  components  at 
each  speed.  For  all  cases  the  helicopter  altitude  was  200  ft  and 
the  sound  level  was  measured  by  a  stationary  observer  located 
near  the  ground  and  at  a  point  250  ft  off  the  flight  path. 

The  flyby  data  referred  to  above  are  compared  in  Figs.  12, 
13,  and  14.  These  figures  contain  data  from  Ref.  4  on  the  meas¬ 
ured  SPL  for  two  runs  at  each  cruise  speed  together  with  a 


/MO 

/» 

/a> 

0-K»III0N 


Run  geometry  for  Figs.  12.  13. 
and  14 
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Fig.  12  -  Sound  pressure  levels  tflyby  data)  for  the  S-5S  helicopter  at  a 
cruising  speed  of  40  knots. 
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curve  calculated  by  the  theory  presented  therein.  In  addition  to 
these  data,  three  curves  computed  by  the  present  theory  are 
included  on  the  figures  corresponding  to  the  first  two  line  com¬ 
ponents  of  each  speed  run.  Two  of  these  curves  represent  calcu¬ 
lated  data  based  on  2.0  and  2.5  inverse  power  laws  lor  the 
harmonic  amplitudes  and  the  assumption  of  random  phase  rela¬ 
tionships  between  the  sound  contributions  of  the  loading  harmonics. 
The  third  curve  computed  by  the  present  method  uses  the  air-load 
amplitude  data  reported  by  Scheiman  for  the  appropriate  cruise 
case,  although  the  phases  of  the  sound  contributions  of  each  air¬ 
load  harmonic  are  again  assumed  to  be  random.  The  helicopter 
data  used  as  input  for  these  calculations  are  given  in  Table  VI 
of  Ref.  4.  \ 

For  the  third  and  fourth  line  components  of  each  speed  run 
only  two  comparative  curves  are  provided  in  the  figures  for  the 
present  computational  method.  The  curve  computed  fros'i  the 
Scheiman  data  is  omitted  because  the  number  of  loading  harmonics 
is  insufficient  to  permit  calculation  of  acoustic  line  components 
above  the  second  for  the  S-58  rotor  system.  Thus,  it  is  to  be 
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Fig.  13  -  Sound  pressure  levels  (flyby  data!  (or  the  S-58  helicopter  at  a 
cruising  speed  of  SO  knots. 


Fig.  14  -  Sound  pressure  levels  (flyby  datai  tor  the  S-58  helicopter  at  a 
cruising  speed  ot  IIP  knots. 


15 


i«nyvuiuMl  jrqrt|  ti  m  V4S 


Kunnc  ot-coMxma  vxm  in 


19.  15  -  Comparison  of  theory  with 
results  c!  whirl-tower  tests  (or  a 
three-bladed  rotor. 


expected  that  the  theoretical  curve  derived  in  Ref.  4  (using 
Scheiman's  data)  would  demonstrate  progressively  poorer  cor¬ 
respondence  with  measured  data  as  higher-frequency  line  com¬ 
ponents  are  considered.  This  obst  vation  is  borne  out  by  the 
cruise-data  comparisons.  For  example,  it  can  be  seen  that  the 
2.0  power  loading-law  assumption  leads  to  significantly  better 
results  for  the  third  and  fourth  acoustic  line  component  data 
(Figs.  12o,  12d,  13c,  13d,  14c,  and  14d).  The  loading-law  com¬ 
putations  include  the  effects  of  harmonic  loading  contributions 
up  to  the  18th  for  the.  third  line  component  and  up  to  the  24th  for 
the  fourth  line  component.  It  is  conceivable  that  the  use  of  this 
many  measured  loading-amplitude  terms  in  the  calculations  would 
lead  to  more  consistently  accurate  results  than  the  use  of  a  load¬ 
ing  ‘Taw.*'  Nevertheless,  it  will  probably  be  some  time  before 
the  blade-loading  characteristics  can  be  measured  or  otherwise 
specified  in  sufficient  detail  to  render  the  use  of  a  loading  law 
unnecessary,  particularly  for  the  calculation  of  the  higher  lii.e 
components  of  the  rotor. 

As  a  final  comparison  of  theory  with  experiment,  calculations 
were  made  to  compare  the  present  rotor-noise  model  with  a  set 
of  results  from  whirl-tower  tests  reported  by  Stuckey  and  Goddard 
(Ref.  13).  The  result  of  this  comparison  is  shown  in  Fig.  15.  Note 
that  calculations  were  made  for  the  case  of  steady  loading  as  well 
as  for  the  cases  corresponding  to  2.5-  and  2.0-power  loading  laws. 
Once  again,  the  result  with  the-  2.0-power  loading  law  is  se  en  to 
be  superior  to  that  with  the  2.5-power  law. 

The  largest  discrepancies  between  theory  and  experiment  in 
Fig.  15  occur  for  the  first  two  acoustic  line  components.  However, 
the  accuracy  of  the  measured  data  at  these  frequencies  is  open 
to  question  and,  in  fact,  Stuckey  and  Goddard  point  out  that  data 
for  the  first,  second,  and  third  line  components  are  suppressed 
due  to  the  limited  frequency  response  of  the  amplifier  and  micro¬ 
phone.  In  addition,  the  authors  also  point  out  that  the  higher  har¬ 
monics  are  exaggerated  due  to  reflections  from  surrounding 
structures  and  to  reingestion  effects.  Thus,  it  appears  that  agree¬ 
ment  between  theory  and  experiment  for  the  case  shown  in  Fig.  15 
may  actually  be  better  than  the  available  data  indicate. 

Note  that  in  the  above  it  has  been  assumed  that  the  effective 
ring  factor  is  0.8  for  every  air-load  harmonic  considered.  To 
be  more  precise,  it  should  be  recognized  that  the  value  of  this 
air-load  factor  is  a  function  of  the  particular  air  load  under  con¬ 
sideration.  However,  since  the  effect  of  this  refinement  on  the 
radiated  noise  is  usually  small  and  data  on  its  specification  is 
normally  unavailable,  it  is  frequently  more  expedient  to  ignore 
it. 
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Fortran  Listing  of  Computer  Program  Implementation 
of  the  Acoustic-Source  Mode! 


The  digital- computer  program  presented  in  this  appendix 
was  written  to  provide  a  rapid  means  of  evaluating  Eqs.  16  and 
19  and  to  produce  punched-card  and  printed  output  specifying  the 
location,  strength,  and  phase  of  the  mcnopoles  comprising  the 
acoustic  array  model  of  ro  or  rotational  noise.  In  addition,  the 
program  provides  for  the  computation  of  the  polar  patterns  (in  a 
vertical  plane)  corresponding  to  the  acoustic  fields  of  each  air¬ 
load  harmonic  contribution  as  well  as  the  total  est;mated  sound 
field  obtained  by  mean-squared  addition  of  these  individual  con¬ 
tributions. 

Specification  of  the  air-load  harmonics  is  accomplished  either 
by  providing  a  nonzero  value  for  the  loading-law  exponent 
(LOADL  7)  or  by  means  of  punched-card  input  that  provide  values 
for  the  harmonic  number  and  normalized  amplitude  and  phase  for 
each  air-load  harmonic  used.  In  the  latter  case  the  parameter 
LOADL. W  is  set  equal  to  zero.  The  highest  air- load  harmonic  to 
be  considered  is  specified  by  the  parameter  NHARM,  and  the 
lowest  (starting)  harmonic  to  be  considered  by  the  parameter 
NSTART.  For  example:  NSTART=1  specifies  consideration  of  all 
air-load  harmonics  beginning  with  the  fundamental  (steady-only 
component)  and  concluding  with  the  harmonic  number  correspond¬ 
ing  to  the  value  of  NHARM. 

Although  the  digital  program  presented  here  is  arranged  to 
provide  polar  pattern  data  for  the  resulting  sound  fields,  it  is  a 
simple  matter  to  modify  it  to  handle  arbitrarily  specified  re¬ 
ceiver  points.  In  fact,  such  a  modification  was  used  in  the 
calculation  of  the  flyby  cruise  data.  In  this  case  data  cards 
giving  the  cartesian  coordinates  of  the  receiver  points  (XR,  YR, 
and  ZR)  were  used  instead  of  the  algorithm  contained  between 
statement  number  19  and  number  21. 

If  it  is  desired  to  account  for  the  phases  of  the  harmonic 
loading  terms,  it  is  necessary  to  sum  their  contributions  to  the 
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sound  field  directly  rather  than  in  the  mean-squared  sense.  The 
printed  output  provides  the  magnitude  and  phase  of  each  air-load 
contribution  so  that  this  computation  may  be  carried  out. 

At  the  conclusion  of  the  program,  control  will  return  to  the 
first  READ  statement,  and  a  new  set  of  input  data  cards  will  be 
processed.  Thus,  to  terminate  the  program  execution,  a  blank 
card  should  be  inserted  after  the  last  set  of  data  cards. 
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THIS  PROGRAM  COMPOTES  THE  LOCATION'  STRENGTH  AND  PHASE  OF  THE 
ACUOSTiC  POINT  SOORCE  APPROXIMATION  TO  THE  FORCF  DIPOLE  FIELD 
OSING  the  EFFECTIVE  RING  PROPELLER  MODEL 
J.MACALOSO  ORL  VERSION  5  9/17/70  MOD  03 

PRINCIPAL  PROGRAM  SYMBOLS  AND  UNITS  ARE  - 

QUANTITY  SYMBOL  ONI TS/REMARKS 


MAGNITUDE  OF  BLADE¬ 
LOADING  HARMONIC  AriARM(l) 

BLADE  ANGLE  AT  K— R  ALPHA 

PHASE  anglf  OF  blade¬ 
loading  HARMONIC  APHSEII  ) 

fld  pt  vert  angle  beta 

CHORD  AT  KR  C 

vflocity  of  SOUND  CO 

CONVERSION  factor  cvfctr 

ANGULAR  SPACING  OF 

FIFLU  points  dbfta 

MONOPOLE  SPACING  DELS 

FREQ  OF  M'TH  HARM 

LINE  COMPONENT  FNDERQ 

RADIAL  FORCE  COMPONENT  FRAD 

MAG  OF  TOTAL  BLADE  FORCE 
AT  K-R  FRCMAG 

DISK  ROTATION  RATE  EREQO 

fld  point  hori/  angl-  gamma 

ROTOR  HEIGHT  H 

EFFECTIVE  KING  FACTOR  < 

MAG  OF  LOAD  LAW  EXP  LOADLW 


ORDER  OF  HARMONIC  M 

HUH  CONVECTION  MACH  NMHK  PACH 

NUMBER  OF  BLADES  NB 

no  OF  FIFLQ  POINTS  NFP 

NUMBER  UE  BLADE 

airload  HARMOMRS  nharm 

FLAG-POLAR  plot  data  nplot 

FLAG-SOURCE ( I ,J)  PNCHFO  NPNCH 

FLAG-SOURCE ( I <J)  PRINT  NPRNT 

MO. SOURCE  DIPOLES  NSOPLS 

LOWEST  LOADING  HARM.  NS TART 

SHAFT  ANG  VELOCITY  OMFGA 

RADIAL  FORCE  ANGLE  PHIC 

ROTOR  TILT  ANGLE  PSI 

TOTAL  TORQUE  CJ 

(AZIMUTH  FUNDAMENTAL) 


NORMALIZED  BY  STDY  TRM 
RADIANS 

INPUT  IN  DEG 

RADIANS 

FT 

FT/SFC 

PSI  TO  MICROBARS 

INPUT  pv  DEG 
FT 

HZ  (Nl)  DUPPLrk) 

LBS 

LBS 

HZ 

INPUT  IN  l)FG 
ET 

FXD  PT  IN  TENTHS 
USED  WHEN  AJRLO  DATA 
Supplied  (INLY  BY  LOADLW 
OTHERWISE  SET  FO  ZERO 


MOTION  RF  +x  AXIS 


I  F= I  DATA  OUTPUT T  ED 
I  F= 1  DATA  OUTPUTTED 
I  F= 1  DATA  OUTPUTTED 


RAD/SEC 
A  T  AN ( FRAD/T ) 

OFGtPOS.  UP  RE  HORIZ.) 
INPUT  IN  LH-iNCHES 


C  BLADE  RADIOS  K  FT 

C  FIELD  POINT  RADIUS  KFP  FT 

C  RPM  RPM  RFV/MIN 

C  MUNOPOLE  SOURCE  MATRIX  SOIJRCF(I.J)  MIXED  UNITS 

C  TllTAL  THRUST  T  LBS 

C  (AZIMUTH  FUNDAMENTAL) 

C  ROTOR  AZIMUTH  ANGLE  THETA  RADIANS 

C  BLADE  LOAD  REF  ANGLE  ThETAR  INPUT  IN  DEG 

c  HUB  VPLUC I  TY  (CONST  )  V  "■  NOT  S  1  NPL'T  (  SUBSONIC  ) 

C  RECEIVES  COORDINATES  XR,YR,ZR  FT 

C 


o  o 


0! MENS  I  ON  SOURCE (5, 2000) , AHARMt 50 ) , APHSE ( 50 ) , ANSTOT ( 2 ♦ 1  JO ) 

R6AL*8  SOURCE.OIST, DBLP,ARG,AIMAG,REAL,L0A0LW*4,MACH*4,TW0PI 
NAMEL1ST/DATIN/M,NFP,RFP,DBETA,H,T,Q,RPM,NB,K,C0,C,R,L0A DLW, FRAD, 
IV, GAMMA /OATHS E/NSOPLS , PH  I C , MACH 

2, 0, FO, ALPHA, DELS ,DTH£T A, OMEGA, FREOO , FNOFRO , NHARM, NST ART, THETAR.PS I 
C  READ  INPUT  DATA 

1  READI5,2)NB,K,CO,C,R.PSI , THE TAR , LOAOLW, V. NPLOT , NPRNT , NPNCH . NHARM, 
INSTART, NSOPLS 

IF (NB.EO.O)  GO  TO  24 

3  READ! 5 ,4)  M, NFP ,RFP , OBET A , GAMMA, H, T, 0, FRAD, RPM 

2  FORMAT (215»E10.3»6F5. 2,615) 

4  F0RMAT(2I5,4F5.1,4E10.3) 

C  PRINT  CAPTION  AND  INPUT  OATA  FOR  THIS  RUN 

Wk I  TE ( 6 ; 5) 

WRITE(fe.DATlN) 

5  FORMATIlHO, 'COMPUTATION  OF  EFFECTIVE  RING  DATA-VERSI ON  5  917/70 
1M0003*/1H0,' INPUT  OATA  FOLLOWS-'//) 

C  COMPUTE  ANO  IDENTIFY  RUN  CONSTANTS 

P 1=3.1415926 

6  T WO PI=6.2B31R5 30 7179586 D+ 00 
OMEGA  =  (RPM*TWOPI )/60- 
V=V<=1 .689 

MACH=V/CO 
G1=1.-MACH~2 
FRF00=0MEGA/TW0P I 
FNOFRO=  FRE00-M5NB 
AK  =  K/10. 

wvnmbr=Twopi*fnofro/co 

DBSTA  =  ( T WOP  I / 360 . )  «OBFTA 
PS1 = ( TWOP I / 360 . ) =PS I 
GAMMA=GAMPA  =  P 1 / 1  BO . 

0  =  0/12. 

FO  =  O/IAKSR) 

CVFCTR=  4.45F+05/(2.54=c2) 

FRCMAG=SORT (T«»2+FO»*2+FKAD--2) 

ALPHA=  ATANIFO/T) 

PH1C=ATANIFRA0/T) 

0ELS  =  O.01  =  ( I T  WOP  I  SCO) / ( M-NB-OMEGA ) ) 

NSRCE  =  2=NS0PLS 
DThETA  =  TWOPI/NSDPLS 
C  PRINT  RUN  CONSTANTS 

WRITE(6,DATUSE) 

THETAR  =  THtTARcT WOP  1/360. 

C  I:.  TIALIZF  ANSWER  MATRIX 

DO  3  lO  LL  =  1,NFP 
310  ANSTOT (2, LL)=U. 

IF  (LOAOLW. EO.O.  )  GO  TO  HOH 

C  COMPOTE  MAGNITUDE  OF  LOADING  HARMONICS  IF  LO.OING  I.AW  IS  USED 
IF(NSTART.NF.l)  GO  TO  1102 
AHARM ( 1 ) = 1 • 

APHSEI 1 ) =0 . 

NSTRT=2 

WR  I  T  E ( 6 , 1 1 50 ) 

1150  FORMAT!'  LAKB0A=1« , 5X , • AHARM= 1 • , 5X , * APHSE=0. • ) 

GO  TO  1101 
1102  NSTRT=NSTART 
1101  00  llOO  NN=NSTRT , NHARM 

AHARM ( NN)  =!./(( NN  -  1 ) ==LOADLW) 


1100 

1151 


APHSE!NN)=0. 

WR1T£!6,:151)NN,AHARM(NN) 

FORMAT ( 1  LAMPDA= • , I 5,5X , 1 AHAKM= • .E10.3,5X,,APHSE=0.' ) 

GO  TO  1109 

READ  £  PRINTOUT  RECORD  OF  AIRLO  HAP.M  DATA-NOTE  NMBR=  1  IS  STEADY 
LOAOING  COMPONENT 
1108  DO  1 2 A  NN=NSTART ,NHARM 

REA0!5,108*  NMBR  *  AH ARM I NN) , APHSE ( NN) 

FORMAT(15,5X,2E10.3> 

WRI TE(6, 150)NMBK,AHAPM(NN) ,APHSF(NNJ 
FORMAT ( 1H  , ' NMBR  =  HARM+1= ■  , I 5,  5X, 1 AHARM=  * , E 10. 

1'  } 


108 


150 


3, • APHSE=  * . E10.3, • DEG 


124 


1 10-* 


APHSE  (  NN)=  APHSE  (NN)CTWiTP  1/360. 

COMPUTE  EO'JI  VALENT  MONOPOLF  SOURCE  DISTRIBUTION  FOR  H0R12  ROTOR 
DO  300  IHM=NSTART .NHARM 
NNM1 = I HM— 1 
DO  8  J=1 .NSRCE .2 
MULT  =  1 
KK  =  J 

THETA  =  -DTHETA  *  i (KK+1 ) /2) *DTHETA 
XXX  =  APHSE  t I HM) 

FMAG=FRCMAG*AHARM( 1 HM ) s ( COS  1 NNM1 *  I THET A- THET AR ) -XXX ) ) 

SOURCE ( 1 ,KK  )  =  AK*R*COS( TnETA)+MULT*Sl N( ALPHA ) »SI N( THETA )* 

1  ( DELS/2. J-MULT-COS ( THETA 1»C0S( ALPHA) -SI NI PHI  Cl- (DELS/2. ) 

SOURCE ( 2 .KK )  =  H+ MULT*COS< ALPHA) *<DELS/;.>*C0S(PHIC> 

SOURCE! 3, KK)  =  -AK-R-SI N( THE fA) +MULTSSI NIALPHAI-COS ( THETA )* 

1  (OELS/2.)  +MULT*S I N(  THETA )*C«)S(  ALPHA )*SI  N(  PHI C  )  =»  I  DELS/ 2.  ) 
S0URCF<4,KK)=MIIlT*SRCMAG( fkag,thopi .DELS. DTHETA) 

SOURCE ( 5 .KK )  =  PHASEITHFTA, OMEGA, AK.R.C.M.NB) 

IF(MULT.LT.O)  GO  TO  8 

MULT  =  -1 

KK  =  KK  +  1 

GO  TO  9 

CONTINUE 

IFIPS1.E0.0.)  GO  TO  330 

DISTRIBUTION  FOR  tilt  ANGLE  IF  ROTOR  NOT  HOR 1 2 


331 

330 
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CORRECT  SOURCE 
DO  331  J= 1 , NSRCE 
X2Ru=bUURCE ( 1 , J ) 

Y2RO=SOURCF(2,J) 

SOURCE! 1,J)  =  XZR0*C0S(PSI ) - ( Y2R0-H) * S I N( PS  I ) 

SOURCE  I  2,  J)=X2R0-S1N(PSI  )  +  (  Y2R0-HJ<=cnSiPSI  )+H 
I F (NPNCH.NE. 1 )  GO  TO  13 

PUNCH  EOIV  ACOUSTIC  SOURCE  OATA  CDS  FOR  CURRENT  AIRLO  HARM 
DO  11  J=l, NSP.CE 

(SOURCE! I ,J) .1=1 .5) 


12 

11 


.R) 


212 
i  ^ 


WRITE (7, 12) 
FORMAT ( 5E15 
CONTINUE 
WRI TE(7,?12) 
FORMAT (ROX ) 

I  F  !  NPP.NT  .NE  . 
PRINT  EOUIV 
WRITFI6.16) 

00  15  J=l, NSRCE 


i)  c-o  TO 
ACOUSTIC 


18 

SOURCE 


DATA  AS  OUTPUT 


FORMAT { 1H0, ' SOURCE  CARDS  DEVELOP  FO  FnR  FFFtCTIVE  KING  APPROX  TO  PR 
10PFLLER  ARE- ' / 1H0 ,,CAR^'.T16,,XI^FT),»T36.,YI!F^),.TS6»,2I(FT),. 
2T71 , • MAGNITUDE ( LB /SOFT ) • . T9 3» • PHASE ( RADI ANS )  •  ) 

17  R OR MAT ( 1h  , I  5,5! 5X, El  5.8) ) 

15  CONTINUE 

18  IF(NPLOT.NE.l)  GOTO  24 


21 


c 


s.^*. 


PRINT  DATA  FOR  POLAR  PLOT  OF  ACOUSTIC  FLO  OF  CURRENT  AIRLD  HARM 
WRI TF ( 6  « 1 19)  1 HM 

19  00  20  L  =  1 » NFP 
RFTA  =  (L-l)’OBFTA 

xr  =  rfp«cos(B6ta)*cosigamha; 

YK  =  H-RFPsSIN(RETA) 

ZR=RFP-COS(BETA)sSI N (GAMMA ) 

RFAL  =  0. 

A 1  MAG  =  0. 

21  00  22  N= 1 ,NSRCE 
XIJPT  =  SOURCE  (  1  ,N)-XR 

RSQRf)=  1 VR— SOURCE  1 2  »  N)  )  «=s=2+  1  ZR-SOURCEI  3,  N)  )s*2 
01 ST= ( -MACHSXOP  TtSORT ( XOPTss^+GlsRSORO) ) /G 1 
0 1 S  T2=  D1 S  T«-2 
XF  =  X0PT-MACH»()1  ST 
F2=ABSI 1.-RS0R0/0IST2) 

351  CSANGL=  S1GN(S<JKT(F2)  .XR-XE) 

ORLP=  (M«NR3l|MfcGA)/C0 
AKG= SOURCE (5.N)-0HLP»D! ST 
ARG  =  l)MOO  (  ARG «  T  WOP  1  ) 

F1  =  1.-MACH=»CSANGL 

COKCTNs  1  .+  I  (  (CSANGL-MACH)SMACH)  /  (  WVNMRK-OI  ST- F 1  )  )«2 
PS  =  0  . 

I F ( CORCTN.GT .0 . )  GO  TO  350 
PS=P!/2. 

CORCTN=-CORCTN 

350  CORCT  N=SOR T ( CORCT N) / ( F 1--2 ) 

OISTl=OIST 

COKANG  =  AT  AN ( ( { MACH-CSANGL ) -MACH ) / ( WVNMt)R»DI SI 1-F1 )  ) 
CORANG=AMOO(CORANG,2.*PI ) 
arg=arg+COkamo,  +pS 

REAL  =  REAL+  ( SOURCE ( A«  N) sOCOS ( ARG )-CORCTN)/OIST 

22  A1MAG  =  AIMAG+  (  SOURCE  (  4»fJ)  —  OSI  N  (ARG )  ^CIiRCTN)  /  OI  ST 
SOMMAG  =  OSUK  T (KFAL--2  +  AI MAG --2 ) / ( SOR  T ( 2 . ) ) 

IF (AIMAG.NF.O.)  GO  TO  222 
I!  (RFAL.NF.O. )  GO  TO  222 
SMPHSF=0. 

G(.  TO  22  3 

222  SMPHSF=  0ATAN2 < A  1  MAG ,KFAL ) 

223  OEGBFT  =  BF T A - ( 360. /T WOP  I) 

DEGPS  =  SMPHSF=( 360./TW0P1  ) 

AMPPSI=SOMMAG/lA4. 

AS0R0=AMPPSI--2 
I F { AMpPSI .NE.O. )  GO  TO  323 
AMPOB=  -1.0F+10 
GO  TO  23 

323  AMPOH=  20.=AL0G10(AMPPSI*CVFCTR) 

23  WKI TF(6,122)  L,R FP, OEGBFT, SOMMAG, SMPHSE,OFGPS, AMPPSI .AMPOB.ASORO 
122  FORMAT (1H  , I5,T11,F10.3,T26,E10.3,TA1,E10.3,T56,E10.3,T71,E10,3, 

1Tm7,F1C.3,T101,c10.3,T116,F10.3) 

ANSTOT ( 1 ,L )=OEGBET 
ANSTOT (2,L)=ANSTOT(2,L)+ASORO 

20  CONTINUE 
300  CONTINUE 

C  PRINTOUT  POLAR  PATTERN  OATA  OIlF  TO  FFFFCTS  OF  ALL  AIRLO  HARM  OSEO 

119  FORMAT ( 1H0, ' POLAR  PATTERN  DATA  FOR  LAMBDA  =  * , I  3 , *  FOLLOWS' 

1/1H  ,'FLO  PT' , Til , 'RADIUS  1  F T ) • , 

1  T26,' ANGLE (DEG) ' ,TA1, ' AMPLITUDE' , T 56, • PHASE ( RAO) ' ,T7l , 'PHASE 

21  DEG) ' ,TR6, 'AMP LIPS! ) ' ,T103,*A(DB  1 UB ) ' , T 1 1 fl» ' A**2  PS!') 

WRI TE(6,311 ) 


-S'  vj^d.  - . .  <->  :-  . 


311  FORMAT 11HO.’ POLAR  PATTERN  DATA  FOR  TOTAL  PRESS  DUE  TO  ALL  AIKLO  HA 
1RM  FOLLOWS- 1 / 1H  ANGLE ( OEG) 1 »T2 1. *  TOT  PRESS! PS  1) 1 , T41 ,• TOT  PRESS! 
20B  RE  1  OB) * ) 

DO  312  LL= 1 « NFP 

ANSTOT(2,LL)=SORTCANSTOT!2,LL) i 
ANSOB=20.’ALOG10!ANSTOT(2,LL)*CVFCTR) 

312  WRITE<6,313)ANST<>T(1,LL) , ANSTOT ! 2 . LL ) , ANSDB 

313  FORMAT ( 1H  ,  E 10 . 3 , T2 5, ElO. 3, T45, E 10. 3) 

GO  TO  1 

24  WR1 TE ! 6 »2 5 ) 

25  FORMAT (1H0,' PROGRAM  COMPL FT ED-NORMAL  EXIT') 

STOP 

END 

THIS  FUNCTION  SUBPROGRAM  COMPUTES  THE  PHASE  OF  ANY  SOURCE  LOCATED 
ON  THE  PERIPHERY  OF  THE  EFFECTIVE  RING  OF  THE  ROTOR. 

REAL  FUNCTION  PHASE’S . THETA , OMEGA, AK, R , C , M, NB ) 

PHASE  =  -  (  (  THETA/ OMEGA  )  +  ( C  /  1  2.  ’OMEGA’AK’R  )  )  )  $M’NB’ OMEGA 

RETURN 

END 

THIS  FUNCTION  SUBPROGRAM  COMPUTES  THE  MAGNITUDE  OF  THE  FOURIER 
COMPONENTS  OF  THE  SERIES  EXPANSION  OF  THE  FORCE-TIME  CURVE 
RELATING  TO  THE  EFFECT. VE  RING  APPROX  TO  THE  ROTOR  NOISE  FIELD 
THIS  VERSION  ASSUMES  AN  IMPULSIVE  PRESSURE-TIME  HISTORY 
REAL  FUNCTION  SRC MAG’8! FRCMAG.T WOPI. DELS. DThETA) 

REAL’S  TW0P1 

SRCMAG=  I FRCMAG/ ( TWOP I  ”2*0ELS) J ’DTHETA 

RETURN 

END 

’  THIS  IS  A  SLASH  ASTERISK  CARD 


f/L« WUrti4li\t|. tisfj AJ  K^fULM.iiVy  tf *>\  Vi .  i  lltf  til VYiVMJ AM *> 
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